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We report inelastic and elastic neutron scattering, magnetic susceptibility, and heat capacity 
measurements of polycrystalline sodium ruthenate (Na3Ru04). Previous work suggests this material 
consists of isolated tetramers of 5* = 3/2 Ru^^ ions in a so-called "lozenge" configuration. Using 
a Heisenberg antiferromagnet Hamiltonian, we analytically determine the energy eigenstates for 
general spin S. From this model, the neutron scattering cross-sections for excitations associated 
with spin-3/2 spin-tetramer configurations is determined. Comparison of magnetic susceptibility 
and inelastic neutron scattering results shows that the proposed "lozenge" model is not distinctly 
supported, but provides evidence that the system may be better described as a pair of non-interacting 
inequivalent dimers, i.e double dimers. However, the existence of long-range magnetic order below 
Tc ~ 28 K immediately questions such a description. Although no evidence of the lozenge model is 
observed, future studies on single crystals may further clarify the appropriate magnetic Hamiltonian. 



I. INTRODUCTION 

Magnetic materials have received continuous research 
interest since the initial description of the Heisenberg 
Haniiltonianiii^i^iii. This is due to interest in both the 
possible technological impact^, as vsrell as, fundamental 
physical phenomenon that many such materials displaj^. 
Pcrovskite-based alkali metal ruthcnates have just re- 
cently started gaining aXteiAiov^^i^^. The ruthen- 
ates exhibit a range of properties from ferro- and para- 
magnetism to superconductivityiiii^iii^, and have been 
shown to demonstrate an interesting cross-road in con- 
densed matter physicsi^ii^ii^iii. The Na-Ru-0 system 
has sparked interest into the many different analogs 
which present various magnetic properties from short 
magnetic order to paramagnetic behavior.— Na3Ru04 is 
one such analog that has induced similar queries over it's 
magnetic structure. 

The structure of Na3Ru04 was first examined by Dar- 
riet et alt^. The refinement of the crystal structure 
in these preliminary measurements show that Na3Ru04 
consists of oxygen coordinated sodium and ruthenium 
sites within the ab plane, separated by a single layer 
of sodium sites displaced along the c-axis. The struc- 
ture of Na3Ru04 was recently re-refined and was deter- 
mined to be mono-clinic, with space group C2/m and 
lattice parameters a = 11.0295(6) A, 6 = 12.8205(7) A, 
c = 5.7028(3) A, and /3 109.90(3)°iaii^. Figure [T^a) il- 
lustrates a single plane of Ru ions together with coplanar 
oxygen and sodium ions. The Ru ions are octahedrally 
coordinated through shared oxygens, and each Ru^^ ion 
may be modeled as having a local spin S* = |. This 
arrangement of ions suggest a local tetramer or lozenge 
spin system as shown in Fig. [Ha)-(b). An isolated spin- 
lozenge with exchange constants J = 3.36 meV and 



a J = 3.88 meV was first proposed by Drillon et al. in 
order to describe magnetic susceptibility measurements 
on this material's. However, the existence of long range 
antiferromagnetic order below T » 30 K was estab- 
lished using Mossbauer spectroscopy^!, and provided the 
first indication that the suggestion of antiferromagnetic 
tetramer clusters in Na3Ru04 may be incorrect. Recent 
measurements of magnetic susceptibility have also been 
interpreted in terms of a spin-tetramer modeliS. Tem- 
perature dependent neutron diffraction studies have con- 
firmed the existence of long-range magnetic order below 
T « 30 K. This long-range order immediately calls into 
question the accuracy of an isolated spin tetramer model 
of Na3Ru04. 

In the following sections, we present the exact ana- 
lytical solutions for the energy eigenstates of a general 
Heisenberg spin S tetramer, and then apply these cal- 
culations to the case of spin-3/2 to determine the exact 
zero-field magnetic susceptibility and inelastic neutron 
scattering intensities. Through an examination of ther- 
modynamic and spectroscopic properties of Na3Ru04 
and a comparison with theoretical predictions for vari- 
ous isolated spin tetramers, we shall see that the isolated 
tetramer model is indeed inappropriate for Na3Ru04. 



II. SPIN-S COUPLED DIMER MODEL 

Figure [ijb)-(d) shows the individual configurations of 
the coupled dimer models we examined. As discussed 
above, Na3Ru04 is suggested to consist of isolated Ru^+ 
(spin-3/2) tetramer clusters, where the tetramers arc in 
a lozenge configuration. This is illustrated in Fig. [TJb). 
In this case, four Ru^+ ions have super-exchange inter- 
actions through Ru-O-Ru bonds, where the bond lengths 
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FIG. 1: Crystal structure and potential exchange interactions 
for Na3Ru04. (a) Crystal structure as viewed along the c-axis 
showing a single plane of atomic sites for ±| unit cells. Ru, 
O and Na sites are black, gray, and white respectively. The 
monoclinic C 2/m unit cell is shown as dashed grey lines. 
(b)-(d) are three possible models of exchange interactions be- 
tween Ru sites in the ab plane as discussed in the text. The 
double dimer model is not spatially confined, and represents 
two individual dimers. 



were determined by neutron diffraction and arc given as 
3.20, 3.20 and 5.56 A for distances di, d2, and d^, re- 
spectively, as illustrated in Fig. [3^. A more general case 
of the spin-lozenge model includes a non-zero exchange 
interaction 7 J resulting in the coupled dimer configura- 
tion, Fig.[ljc). If J = 0, one recovers two isolated dimers 
or a double dimer configuration, c.f. Fig. [IJd). We note 
that the double dimer model is not spatially confined to 
the four Ru^+ ions in the lozenge configuration, but could 
also represent other dimer interactions in the Na3Ru04 
crystal structure. Using the coupled dimer model Hamil- 
tonian, we determine the eigenstates for general S*, and 
calculate the corresponding magnetic susceptibility for 
fitting purposes. Then, with the choice of the appropri- 
ate magnetic ground state, the excitations observed with 
inelastic neutron scattering (INS) and their correspond- 
ing structure factors are also determined. 



A. Hamiltonian and Energy Eigenstates 

All three of the dimer configurations in Fig. [T] can be 
described by a single Hamiltonian, Eq. [1] By using this 



model, we can clearly examine the three possible con- 
figurations that may describe Na3Ru04. Using nearest- 
neighbor Heisenberg interactions and a Zeeman magnetic 
field term for magnetic fields B defining the z-axis, the 
general Hamiltonian is 

n = j[(Si • S3 + Si • S4 + S2 • S3 + S2 • S4) 



aSi ■ S2 + 7S3 • S4 



[SI + 5| + 5| + Sl)giiBB, 



(1) 

where aJ is the interaction for the a-dimer, 7 J is the in- 
teraction for the 7-dimer, and /is is the Bohr magneton. 
We define the exchange interaction as positive for anti- 
ferromagnetic interactions, and Si is the quantum spin 
operator for a spin-S" ion at site i=l,2,3,4. The Zeeman 
term interacts with the 2-component of the spin Hamil- 
tonian, lifting the degeneracy of magnetic substates in 
applied magnetic fields. 

The Hamiltonian, Eq. [1] is rotationally invariant in 
spin space, such that the total spin, Stoti and Sz are 
good quantum numbers. For the general case of a spin- 
s' tetramer cluster, the energy eigenstates have the total 
spin decompositions given by 



nS'tot = l 



StotY 



(Stot + l)(Stot+2 



„i(4S+2+8StotS+St, 
"^tot 



(2) 

where the total number of magnetic states in a general 
S tetramer are (25' + l)"*. Therefore for the S = 3/2 
tetramer, the energy eigenstates have S"* = 256 magnetic 
states and the spin decomposition is given by 

|«.|®|®| = 
Individual Spins, 



(3) 



(3 © 2 © 1 ® 0) ® (3 ® 2 © 1 © 0) = 
Dimer States, and 

6 © 53 © 46 © 31" © 2" © 1^ © 0^ 
Tetramer States. 



The superscript in the tetramer states denote multiple 
Stot states. Each multiplet containing 2Stot + 1 magnetic 
states, which are degenerate given an isotropic magnetic 
Hamiltonian such as the Heisenberg form of Eq. [U where 
the degenerate states can be split by a magnetic field. 
This breakdown of the dimer and tetramer spin states 
helps clarify which dimer states are interacting to create 
the composite tetramer states. 

By expanding the Kambe approach2i2^, we can rewrite 
the Hamiltonian in terms of total spin for the individual 
diagonalizablc components, in which the eigenstates and 
eigenvalues of the Hamiltonian may be found by diago- 
nalization in the convenient basis of two dimers. This ap- 
proach gives information about the states of the dimers 
as the tetramer states are determined, which allows a 
clearer picture of the magnetic excitations. In practice. 
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we employ the usual set of z-polarized magnetic basis 
states. The energy levels are then determined simply by 
considering a dimer basis, where Sa corresponds to the 
spin state of the a dimer and Sy corresponds to the spin 
state of a 7 dimer as described in the Hamiltonian, Eq.[T] 
Using this dimer basis, the energy levels for the general S 
coupled dimer can be determined exactly and are given 

by 



E 



J 



[ytot + ^„(a - 1) + ^^(7 - 1) - 2(a + 7),^] (4) 



where .ytot = Stot{Stot + 1) with Stot denoting the mag- 
netic state of the system, S^a = Sa{Sa + 1) and Sa is the 
spin state of the a dimer {S1-S2 dimer), 5^-^ ~ S-y{Sy + 1) 
and S-y is the spin state of the 7 dimer (S3-S4 dimer), and 
J/' = S{S + 1) with S being the spin of the system. 

For the ^ = | tetramer, the energy cigenstatcs of the 
isotropic Heisenberg Hamiltonian with general S^^^ are 
given in Table [J The magnetic substates are all degener- 
ate in the absence of an applied magnetic field, but can be 
split linearly in accordance with the Zccman field term. 
The individual states can be split into 2Stot + 1 states 
with magnetic field. Therefore, the energy levels split by 
gfiBBS^gf, where S^^^ is from . . . (25404-1-1). The ground 
state of the = | tetramer can either have a nonmag- 
netic S' = ground state or a magnetic S = 1 ground 
state depending upon the values of a and 7. Simple ex- 
amination of the energy levels in Table |T] indicates that, 
assuming antiferromagnetic exchange for J, the ground 
state will be non-magnetic when both a and 7 are less 
than |. 



B. Spin- 1 Magnetic Observables 

Magnetic observables, i.e. specific heat and magnetic 
susceptibility, associated with the general coupled dimer 
model can be derived via the partition function and 
eigenvalues. Thus, such macroscopic measurements of 
Na3Ru04 may serve to place limits on the nature of the 
interactions. Wc also determine the excitation energies 
and structure factors which would be observed in INS 
measurements. Both thermodynamic and spectroscopic 
measurements should be consistent for any appropriate 
description of the experimental system. 



1. Magnetic Susceptibility and Heat Capacity 

We now present the method for determining the parti- 
tion function and magnetic susceptibility for the S* = 3/2 
coupled dimer model. Due to the length of the equations, 
we present the eigenstates and eigenvalues explicitly in 
Table [T] and represent the magnetic observables as sum- 
mations over energy eigenvalues. Using this method, the 



TABLE I: Energy levels for a spin-| lozenge 



Stot 



Spin State" 

Stot Stot >Sa,S.y 



I 6 Stot >3,3 



5 Stot >3,3 
5 Stot >3,2 
5 Stot >2,3 



4 Stot >3,3 
4 Stot >3,2 
4 Slot >3,1 
4 Slot >2,3 
4 Sfot >1,3 
4 Sfot >2,2 



3 Stot 
3 Sfot 
3 Sfot 
3 Stot 
3 Sfot 
3 Sfot 
3 Stot 
3 Sfot 
3 Stot 
3 Stot 



>3,3 
>3,2 

>3,1 
>2,3 
>1,3 
>2,2 
>2,1 
>1,2 
>3,0 
>0,3 



2 Stot 
2 Sfot 
2 Sfot 
2 Sfot 
2 Sfot 
2 Stot 
2 Stot 
2 Stot 
2 Stot 
2 Stot 
2 Stot 



>3,3 
>3,2 
>3,1 
>2,3 
>1,3 
>2,2 
>2,1 
>1,2 
>1,1 
>2,0 
>0,2 



1 Sfo 



Sfot 
Stot 
'-'tot 

QZ 

^tot 
Stot 
Stot 

CZ 
^Jt.nt 



1 Sfot 



>3,3 
>3,2 
>2,3 
>2,2 
>2,1 
>1,2 
>1,1 
>1,0 
>0,1 



Sfot >3,3 

Sfot >2,2 

Sfot >1,1 

Sfot >0,0 



Energy 
Level'' 



J(9 + |q+|7) 



J(3 + |a+f7) 



J(6 
J 6 



7) 



ia+ i7) 



j(_2 + |a+|7) 
J(l + |a-|7) 

J(3 + |a-^7) 
J(l-|a+|7) 

J(3-^Q + |7) 
j(4-|a-|7) 



J(-6 + |a+|7) 
J(-3 + |a-|7) 
J(-l-f |a-^7) 
J(-3-|a+|7) 
J(-l-^a + |7) 

J(2-|a-^7) 
I7) 



J(2 
J(fa 



f7) 



h) 



J(-9 + |a+|7) 
J(-6 + |a-|7) 
J(-4-f |a-^7) 

J(_4-fa + |7) 
J(-3-|a-|7) 
J(-l-fa-^7) 
J(-l-^a-|7) 
J(l-fa-ii7) 



7) 



J(-ll + |a+|7) 
J(-8-f |a-|7) 
J(-8-fa+|7) 
J(-5-|a-|7) 
J(-3-|a-f7) 
J(-3-^Q-|7) 
J{-l-4a 



^7) 



J(->-^7) 



J(-12 + |a+|7) 
J(-6-|a-|7) 
J(-2-fa-|7) 



J(-fa-f7) 



""Sa and S-y axe the spin state of the a and 7 dimers. 
''The magnetic field splitting of the states is given by the 
addition of a Stotgi^sB to the energy of the magnetic states. 
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canonical partition function is given by 

N 



(5) 



(6) 



and the magnetic susceptibility is given by 

= Ee. (2Stot + 1) {Stot + 1) Stot e-^'Ei . 

In these formulas, the sum Ei^i over all N indepen- 
dent energy eigenstates (including magnetic substates), 
the sum is over energy levels only, = mgfiB 

where m = S^Qt/fi is the integral or half-integral mag- 
netic quantum number, g is the Lande p-factor, (3 = -^Kj: 

and ks is Boltzmann's constant^. The heat capacity can 
also be determined from 



C = ksP 



2 dMz) 



(7) 



In general, the heat capacity is especially useful for con- 
firming the proper accounting of eigenstates through a 
numerical calculation of the magnetic entropy of the 
tetramer system at large temperature. 



kn 



|21n(2), 
lln(2M), 



dp 

T 



a < I and 7 < 



> I or 7 > 



(8) 



3' 



Here Af is the dimensionality of the full Hilbert space and 
A/q is the degeneracy of the ground state manifold. The 
numerical calculation of entropy within various regions 
of a and 7 confirms the ground state of the spin-3/2 
tetramer. 



2. Inelastic Neutron Scattering 

The experimental focus of this paper is the use of INS 
to investigate the nature of the magnetic interactions and 
their respective excitations. Using methods presented in 
Haraldsen et al^^^, we next determine the excitation 
energies and structure factors for the observable transi- 
tions of the coupled dimer models shown in Figs.[TJb)-(d). 

For transitions out of the ground state, the excitation 
energy, huj, is simply the difference in energy between 
the excited and ground states. Such excitations would be 
non-dispersive in the absence of inter-tetramer exchange. 
The differential cross-section of finite systems is propor- 
tional to the neutron scattering structure factor 



Siq) = \Fm' 



^ {^f{Xf)\Va\^,)' 



(9) 



where F(g) is the magnetic form factor and the vector 
Va{q) is a sum of spin operators over all magnetic ions 
in a unit cell 



(10) 



For rotationally invariant magnetic interactions and an 
Stot — ground state in the T = limit, only Stot = 1 
final states as shown in Tableware observable via INS. 
However, due to the nature of the tetramer states as be- 
ing composite dimer states, this implies that it is only 
possible to excite transitions of the individual dimers 
which make up the tetramer structure, ASq,/^ = ±1,0. 
To interpret neutron experiments on powder samples, we 
require an orientation average of the unpolarized single- 
crystal neutron scattering structure factor. We define 
this powder average by 



(11) 



With respect to the spin-3/2 rhombus model, the val- 
ues of the magnetic interactions quoted in the literature 
suggest a Stot = ground state, with dimer spins 5q = 3 
and S^ = 3 ^°i^° . Therefore, due to this selective restric- 
tion of the spin excitations, only three of the nine Stot = 1 
states are accessible from that ground state through INS. 
The respective excitation energies {Es^^t^ScSj) are 



0,3,3^1,3,3 



= J, 



-^0,3,3^1,3,2 = J(4 - 37), 



(12) 



Eo. 



3,3^1,2 



,3 = J(4-3a), 



and the powder average INS structure factors 
{S{q)stot,Sc,Sj) for these transitions are 



-5(9)0, 



3,3^1,3,3 



2|F(g )|^(2 - 4jo{qdi) + joiqd2) + Mqds)), 



^(9)o,3,3-.i,3,2 = ^^(1-^0(9^3)), 



(13) 



-(1 - k{qd2))-, 



'5'('i')o,3,3^1,2,3 — — 

where di, ^2, and da are the interatomic separations 
(shown in Fig. [T]), jo(x) = ^^^^p^, and |F(g)| is the 
Ru^+ magnetic form factor (a parameterization is given 
by Parkinson et alS^^). The transition of |00 >3,3^ 
\lStot >3,3 is an excitation of the full tetramer, while the 
other two transitions are excitations of individual dimers. 
This shows that out of the nine possible spin-1 states to 
be excited for the spin-| model, only three are accessible 
by INS and will have a corresponding intensity profile. 
The unseen transitions are inaccessible because they re- 
quire multiple spin transitions, and since a neutron can 
only provide one transition, only transitions that excite 
the individual components will be observed by neutron 
scattering. 
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III. EXPERIMENTAL TECHNIQUES 

Powder samples of Na3Ru04 were prepared by solid- 
state reactions from stoichiometric amount of NaOH and 
Ru02. The starting stoichiometric mixture was initially 
ground together and then held at 500 °C for 20 hr. un- 
der an O2 atmosphere. After re-grinding, the powder was 
heated to 650 °C for another 20 hr, again under an O2 at- 
mosphere. The resulting dark grey powder was reground 
and checked for impurity phases using X-ray powder 
diffraction. If any impurity phases were evident, the pow- 
der was refired and the process repeated. This growth 
procedure is similar to that described in Ref. [Tol . Powder 
refinement of room temperature X-ray diffraction mea- 
surements yielded lattice parameters of a = 11.012(7), 
b = 12.809(9), c = 5.687(3) A, and /3 = 109.91(3)° for 
the C2/m monoclinic unit celU^. These values compare 
well to the fully refined structure described in Ref. [lol . 
Single crystals of appropriate mass are unfortunately not 
yet available for INS measurements. 

Heat capacity measurements were performed on a 
small single crystal of mass « 10 mg, which was obtained 
through the synthesis procedure described above. This 
single crystal grew as a small platelet, with the c-axis 
normal to the plane of the platelet. Heat capacity mea- 
surements were performed with a commercial calorimeter 
between T = 1.8 K and 300 K, using the relaxation tech- 
nique. Measurements were carried out in zero and 8 T 
applied magnetic fields, with the field applied along the 
c-axis of the single crystal sample. 

Magnetization measurements were performed on pow- 
der and single crystal samples using a commercial 
SQUID, as a function of applied magnetic field and tem- 
perature. SQUID measurements on the same single crys- 
tal sample that was used for heat capacity measurements 
agree well with those taken on a powder sample. 

INS measurements were performed using the MARI 
time-of-fiight spectrometer at the ISIS neutron scattering 
facilitji^l. The sample consisted of ~ 45 g of Na3Ru04 
powder in a square aluminum foil sachet (approximately 
50 by 50 by 8 mm), suspended from the cold-tip of a 
closed-cycle He'' refrigerator. The sachet was oriented 
with the 50x50 mm surface normal to the incident neu- 
tron beam. An incident energy of Ei = 25 meV was 
used, and data were taken at several temperatures be- 
tween r = 8 K and T = 305 K. This configuration re- 
sulted in a measured instrumental energy resolution at 
the elastic position of Shoj = 0.982(7) meV full width at 
half maximum (FWHM). Data were corrected for detec- 
tor sensitivity through room temperature measurements 
on a vanadium standard. 

INS measurements were also carried out using the HB3 
triple-axis spectrometer at the high flux isotope reac- 
tor (HFIR) at Oak Ridge National Laboratory. For 
these measurements, the sample consisted of 20.7 g of 
Na3Ru04 powder in a cylindrical aluminum sample can 
of 18 mm diameter and 57 mm height. The sample can 
was scaled under He gas and mounted to the cold-tip 
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FIG. 2: Heat capacity of a Na3Ru04 single crystal. Mea- 
surements were performed at zero field (black squares) and 
H — 8 T (red open circles), with H \\ c. Error bars represent 
an estimated five percent error in the measurement. 



of a closed-cycle He^ refrigerator. Horizontal coUima- 
tion was chosen as 48' — 40' — 40' — 120' between source 
and monochromator, monochromator and sample, sam- 
ple and analyzer, and analyzer and detector, respectively. 
The spectrometer was operated with fixed final energy, 
Ef ~ 14.7 meV, using a pyrolytic graphite (PG 002) 
monochromator and analyzer. Pyrolytic graphite fil- 
ters were placed after the sample to substantially reduce 
higher-order spurious scattering processes. In this con- 
figuration, the energy resolution at the elastic position 
was Shoj = 1.10(2) meV FWHM, as measured from the 
incoherent scattering at Q = 1.2 A^^. The wave vec- 
tor resolution was measured to be SQ ~ 0.0407(7) A"'^ 
FWHM using the (110) nuclear Bragg peak. All mea- 
surements were made for fixed incident neutron monitor 
count. 

Elastic neutron scattering measurements were also per- 
formed using the HB3 triple-axis spectrometer, with 
Ei = Ef = 14.7 meV. These measurements were per- 
formed on the same powder sample as the inelastic HB3 
measurements, with horizontal collimation 48' — 20' — 
20' — 70'. This resulted in an energy resolution at the 
elastic position of 6huj « 0.8 meV FWHM. The wave vec- 
tor resolution was measured to be SQ = 0.0254(9) A^^ 
FWHM using the (110) nuclear Bragg peak. 

INS measurements were also performed to place lim- 
its on the value of a possible energy gap in the exci- 
tation spectrum. These were performed using the IRIS 
backscattering spectrometer at the ISIS neutron source 
at the Rutherford Appleton Laborator}*2^. The sample 
measured was the identical powder used for the HB3 mea- 
surements. The IRIS spectrometer was operated at 25 Hz 
with cooled PG002 analyzers (T = 10 K) and a Berylium 
filter (T = 25 K) to avoid contamination from higher- 
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FIG. 3: Magnetic susceptibility of Na3Ru04 powder (black 
squares), showing fits using the three models: Lozenge (solid 
blue), coupled dimer (dashed green), double dimer (dotted 
magneta). The inset shows x~^i with a Curie- Weiss fit as 
described in the text (which gives lOcw-'l = 14.0 meV, and is 
consistent with dominantly antiferromagnetic interactions.) 



order reflections, resulting in a 17.5 /leV FWHM energy 
resolution at the elastic position as measured with a vana- 
dium standard. The high-resolution (back-scattering) 
diffraction banks with a resolution of = 2.5 E-3 were 
also used on IRIS. 



IV. EXPERIMENTAL RESULTS 

Figure [5] shows the heat capacity as a function of tem- 
perature for H = and 8 T. There are two clear lambda- 
like anomalies at T w 23 and T « 28 K, signifying 
phase transitions at these temperatures. Previous neu- 
tron diffraction measurements have shown the existence 
of only a single, broad phase transition near 30 K in 
Na3Ru04, corresponding to the onset of long-range mag- 
netic order—. Our heat capacity measurements indicate 
that the observed broad transition is likely due to two 
successive transitions that occur at similar temperatures. 
No change was noted in these transition temperatures 
when measured at = and H = 8 T. 

The magnetic susceptibility of Na3Ru04 was measured 
over the range 2 < T < 350 K; the resulting data is 
shown in Fig. [3] The susceptibility also shows evidence 
for a phase transition near T « 30 K. The negative in- 
tercept in x~-^(r) and the decrease in x(T) below the 
transition temperature are consistent with dominantly 
antiferromagnetic interactions. 

Figure 3] shows the temperature dependent inelastic 




FIG. 4: INS intensity from Na3Ru04 powder versus energy 
and momentum transfer at (a) T = 8 K, (b) T = 26 K, 
(c) r = 37 K and (d) T = 88 K. These measurements were 
carried out on the MARI spectrometer at an incident energy 
of Ei = 25 meV. 



neutron scattering data taken on the MARI spectrom- 
eter at ISIS. There is significant inelastic scattering in- 
tensity in the vicinity of ftcj w 5 mcV, which decreases 
in intensity rapidly with increasing wave-vector. The 
wave-vector dependence implies that the scattering is 
magnetic in origin. In the T = 8 K data, a weak ex- 
citation near Huj ^ 10 meV is also evident. As the 
temperature increases, the inelastic scattering intensity 
rapidly decreases and moves to smaller wave- vectors, con- 
sistent with an evolution from antiferromagnetic spin- 
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(a) Q = 1 .6 A ' 
□ T = 8 K 
• T = 250 K 
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FIG. 5: INS intensity from Na3Ru04 powder versus energy 
transfer at (a) Q = 1.6 A"^ at T = 8 and 250 K and (b) 
integrated between 0.4 A"^ < Q < 1.6 A"^ at T = 8 K, 30 K 
and 88 K. Black lines are Lorentzian fits, as described in the 
text. The data in (a) are from HB3, and the data in (b) are 
from MARL 
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FIG. 6: (a) INS intensity versus temperature at hui — 18 meV 
and Q = 1.6 A~^. (b) INS intensity versus momentum trans- 
fer at hu = 18 meV at T = 8 and T = 250 K. Data were 
acquired using IIB3. 



waves to paramagnetic scattering with increasing tem- 
perature. We speculate that the excitations observed be- 
low Tjv arc acoustic and optical spin-waves associated 
with the long-range ordered phase. Higher incident en- 
ergy measurements were also performed, which show ev- 
idence for phonon excitations above 20 mcV. 

In Fig. [5l we show the scattering intensity as a function 
of energy transfer for the single wave- vector Q = 1.6 A^^ , 
measured on the HB3 spectrometer, as well as the inte- 
grated scattering intensity for 0.4 < Q < 1.7 
on the MARI spectrometer. Single Lorentzian fits to the 
low-temperature data suggest modes at tku = 5.03 ±0.08, 
9.8 ± 0.2 and 17.9 ± 0.3 meV (for the data shown in 
Fig. IHa)), and hw = 4.95 ± 0.04 and 9.8 ± 0.1 meV for 
the data shown in Fig. [IJb). The increase in intensity 
of the 18 meV excitation with increasing temperatures 
suggests that it is likely a phonon excitation. 

The 18 meV phonon excitation was investigated more 
carefully as a function of temperature and wave- vector, as 



shown in Fig. [SI These data show a monotonic increase in 
scattering intensity as the temperature is increased, and 
an increase in scattering intensity with increasing wave- 
vector. These results further support the identification 
of the 18 mcV mode with a phonon excitation. 

We also examined the temperature dependence of the 
elastic scattering in the vicinity of Q ~ 1 A~^. Fig- 
ure [7] shows the scattering intensity observed in the HB3 
Na3Ru04 powder measurement as a function of tem- 
perature and wave-vector. As the temperature is de- 
creased, there is an increase in the scattering intensity 
at Q ~ 0.99 and w 1.07 A~^ corresponding to the tran- 
sition to long range magnetic ordering. Below T « 25 K, 
the magnetic Bragg peaks appear to move as a func- 
tion of decreasing temperature. This is also evident in 
Fig. [SI which shows the scattering intensity as a function 
of temperature for various individual momentum trans- 
fers. For certain Q values, the scattering intensity shows 
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FIG. 7: Elastic scattering intensity of Na3Ru04 as function 
of momentum transfer and temperature. Data were acquired 
using the HB3 spectrometer. The contour hues correspond to 
intervals of 20 counts per ten seconds, and are only plotted 
for count rates between 80 and 200 counts per ten seconds. 
The data were obtained at 20 temperatures between T = 8.5 
and 34 K. 



non-monotonic temperature dependence, for example as 
shown in Fig. [7] at Q = 0.99 A^^ and 1.05 A'^ These 
behaviors may be due to the presence of two magnetic 
phase transitions near 25-30 K as seen in the heat capac- 
ity measurements. 

High resolution backscattering measurements investi- 
gated the magnetic spectrum for energy transfers below 
1.7 mcV (shown in Fig. [5][a)). Excitations arc populated 
out of the Q ~ 1.1 A~^ wave- vector with a very steep dis- 
persion at such low energy transfers. Figure, ^h) shows 
the difference of the elastic scattering intensity between 
the disordered and ordered phases as measured using 
IRIS and HB3 illustrating the low-temperature powder 
magnetic Bragg peaks in Na3Ru04. These data illus- 
trate that the excitations arc dispersing directly out of 
the magnetic Bragg peaks. The data in Figs. [D^a) show 
no indication of a gap in the magnetic spectrum down to 
- 250 ficY. 



V. DISCUSSION 

The high temperature magnetic susceptibility fitted 
using the finite cluster models for T > 30 K (where there 
is no long range order). Table HIl shows the magnetic in- 
teractions determined from these fits for the three mod- 
els considered. We also compare the extracted magnetic 
interaction parameters to the Curie- Weiss temperature. 
The Curie- Weiss temperature in a mean field approxima- 
tion is 

_ S{S + 1) 
'ow — 2 0' ^ ^ 




T(K) 

FIG. 8: Elastic scattering intensity from Na3Ru04 powder 
versus temperature for several wave-vectors. The data corre- 
spond to the individual wave-vectors shown in Fig. [7] The 
data have been offset along the vertical axis for presentation. 



where, in this case, S = | and Jq is the sum of the mag- 
netic exchange constants^. As shown in Fig. [3l all three 
cluster models qualitatively reproduce the high temper- 
ature susceptibility data, and only deviate strongly close 
to the transition temperature. However, a comparison 
of the calculated Curie- Weiss temperatures based upon 
the double dimcr model is more consistent with the ex- 
perimental Curie- Weiss temperature, \Qw\ — 14.0 nieV 
(antiferromagnetic), illustrated in the inset of Fig.[3l Be- 
cause the Curie- Weiss temperature is proportional to a 
sum of exchange constants, the presence of inter-cluster 
exchange could significantly effect the estimate value of 
&w- For example, there are 28 bonds with distances of 
five to six angstroms between the Ru sites in one cluster 
and the Ru sites in all neighboring clusters. If inter- 
cluster interactions are large, the estimated Curie- Weiss 
temperature would deviate significantly from the values 
quotes in Table HIl 

The inelastic neutron scattering energies and intensi- 
ties can also be calculated using the exchange values de- 
termined from the magnetic susceptibility f Table HH) and 
the ion positions given by Regan et ali^ The three mod- 
els predict INS observable energy gaps of 2.96, 3.08, and 
12.16 meV for the lozenge model, 1.33, 3.49, and 23.53 
meV for the coupled dimer model, and 2.47 and 9.16 meV 
for the double dimer model. We have already noted that 
the data shows magnetic excitations at approximately 
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FIG. 9: (a) T = 15 K Powder INS scattering intensity from 
Na3Ru04 as a function of energy and wave- vector transfer 
measured using the IRIS spectrometer at ISIS, (b) Intensity 
as a function of wave- vector transfer for data acquired at both 
IRIS and HB3. IRIS data were acquired at T = 8 and 30 K 
using the high resolution diffraction banks. HB3 data were 
acquired in same configuration as Fig. [7]at T = 8 and 30 K. 
Both data sets indicate clear peaks at 1.0 and 1.07 A~^. Panel 
(c) shows the high resolution IRIS data fitted to a series of sin- 
gle width Gaussian peaks at g = 1.006(1), 1.025(4), 1.065(1) 
and 1.078(1) A'^ 



5.0 and 9.8 meV, c/. Fig. [5l Figures [4] and [9] also show 
evidence for a spin- wave emerging from the antiferromag- 
netic Bragg peak in the vicinity of the Q « 1 for 
T < Tn. Although a clear transition to long-range mag- 
netic order is evident at low temperatures, it is nonethe- 
less reasonable to examine the properties of these exci- 
tations as observed in INS especially their wave-vector 
dependence, since these are characteristic indicators of 
the nature of the interactions. 

On comparing the predicted excitation spectrum, us- 
ing the magnetic susceptibility data to determine the 
interaction strengths, the double dimcr model seems to 



give the most realistic description of the excitation ener- 
gies. We note however that the inter-cluster interactions 
which produce long-range order may significantly affect 
the energy levels. All three models predict an INS visi- 
ble excitations below 3 meV and above 9 meV. We also 
examine the wave- vector dependence of these three mod- 
els; Figure [TU] shows constant energy scans performed 
above and below T/v at 5.0 and 9.8 mcV energy transfer. 
For comparison, we first calculate the wave-vector de- 
pendence of the INS scattering intensity using Eq. [T3] for 
the lozenge model geometry. These lineshapes (shown in 
Fig. nop are unable to account for the initial rapid rise in 
scattering intensity at small wave-vectors, which implies 
that there are significant exchange interactions between 
spins at larger separations than are present in the lozenge 
model. If the ionic distances are allowed to vary freely 
in fitting the data, an interesting result emerges. For 
the 9.8 mcV excitation, the fitted dimer separation in a 
single dimer model is 5.66±0.10 A bond. This separa- 
tion agrees with the length of the 7-dimer. This makes 
it unlikely that the coupled dimer model is realistic, as it 
predicts a 23 meV 7-dimer excitation using the magnetic 
susceptibility parameter. The lozenge and double-dimer 
models with susceptibility-fitted parameters predict 7- 
dimer excitations at 12 meV and 9 meV, respectively. 
The double dimer model is evidently closer to the ob- 
served gap of 9.8 meV. Fitting the 5.0 meV data with a 
free dimer length in which, both the lozenge and dimer 
models gives a length of 4.60±0.02 A, which does not 
correspond to any Ru-Ru separation in the structure of 
Na3Ru04. We conclude that the intercluster interactions 
are important enough to modify the energies and wave- 
vector dependences of the excitations, so that the three 
simple dimer models do not give a decent description of 
the excitations. 



In Table [TTl gives estimated magnetic interaction pa- 
rameters for the various models using the observed neu- 
tron scattering excitations at huj w 5.03 and 9.8 meV as 
input. Since no third magnetic excitation was observed 
to 30 meV, the coupled dimer model cannot be uniquely 
constrained. Although the lozenge model also predicts 
three excitations, two of those excitations involve dimers. 
Therefore, we can determine both exchange constants 
using the neutron scattering results. Using these ex- 
change interactions, the lozenge model gives an estimated 
Curie- Weiss temperature of 29.5 meV, more than dou- 
ble the observed value from magnetic susceptibility. The 
double dimcr model gives a value of |8iu| = 18.5 meV; 
this is closer to the measured value, although it is still 
30% larger than the observed value (from the magnetic 
susceptibility). If the excitations are indeed well de- 
scribed by the double dimer model, then they demon- 
strate that the 5.0 and 9.8 meV modes correspond to ex- 
citations of the a- and 7-dimers, with the nature of the 
a-dimer excitation being significantly modified coupled- 
cluster effects and the onset of long-range magnetic order. 
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TABLE II: Exchange interactions estimated from the 
magnetic susceptibility and inelastic neutron scatter- 
ing. 



ticipatc that neutron diffraction measurements on single 



Curio- Weiss Fit 
|e„| = 14.0±0.2 meV 


Magnetic Susceptibility 


Model 


J (meV) aJ (meV) 7J (meV) 


|e„| (meV) 


Lozenge 


3.19±0.01 3.11±0.03 


19.8±0.1 


Coupled 
Dimer 


3.7±0.1 4.6±0.3 -4.8±1.2 


18.3±1.3 


Double 
Dimer 


2.69±0.01 9.75±0.05 


15.6±0.1 


Inelastic Neutron Scattering" 


Model 


J (meV) aJ (meV) 7J (meV) 


|e„| (meV) 


Lozenge 


5.03±0.08 3.5±0.1 


29.5±0.1 


Double 
Dimer 


5.03±0.08 9.8±0.2 


18.5±0.1 



"There is insufhcient information to determine the 
interaction strengths in the coupled dimer model. 



VI. CONCLUSIONS 

In summary, we have given analytical results for the en- 
ergy eigenvalues and eigenstates for coupled dimcrs with 
general spin S ions. We derive analytical closed-form 
results for magnetic susceptibility and inelastic neutron 
scattering excitation functions and their wave-vector de- 
pendences for several tetramer models, and compare our 
results to experimental data on the S = 3/2 tetramer 
spin lozenge candidate, Na3Ru04. 

On considering the observed magnetic susceptibility 
and inelastic neutron scattering data and comparing 
these results to several tetramer models, we first find 
that the Na3Ru04 data is not consistent with the spin- 
lozenge model of Regan et alt^. Although no isolated 
tetramer model is able to describe all of the thermody- 
namic and spectroscopic measurements simultaneously, a 
double dimer model, with bond lengths of 4.60 A and 
5.66 A, does provide a description of some aspects of 
the observed thermodynamic properties and the inelas- 
tic neutron scattering measurements on the two observed 
magnetic excitations. However, only one spatial distance 
corresponds to length seen in the material shown in Fig. 

m 

Measurements of the heat capacity and clastic neutron 
scattering data show that there are two distinct magnetic 
phase transitions in this material, at T w 23 and 28 K 
(This is the first evidence for two low-temperature phase 
transitions in this material). Clearly, an understanding 
the nature of these long-range ordered magnetic phases 
will provide useful additional information regarding the 
nature of magnetic interactions in Na3Ru04. We an- 




FIG. 10: a) INS intensity versus momentum transfer of the 
9.8 meV excitation in Na3Ru04 at T = 8 K(blue squares) 
and 30 K (red circles); the black solid lines are fits to the 
dimer model and the gray dashed lines are predictions given 
fixed physical distances. b)INS intensity versus momentum 
transfer at T = 8 (blue squares) and 30 K (red circles) for 
the 5.0 meV excitation. The solid black line is the lozenge 
model fit, the dashed gray line is the dimer model fit, and 
dotted light gray line is a prediction using physical distances 
(The data were taken using the MARI spectrometer at ISIS, 
as described in the text). 



crystal samples of Na3Ru04 will be the most useful next 
step in the experimental studies of this material. 
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